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Abstract. Quantum creation of the universe is described by the density matrix 
defined by the Euclidean path integral. This yields an ensemble of universes - 
a cosmological landscape — in a mixed quasi-thermal state which is shown to be 
dynamically more preferable than the pure quantum state of the Hartle-Hawking type. 
The latter is suppressed by the infinitely large positive action of its instanton, generated 
by the conformal anomaly of quantum matter. The Hartle-Hawking instantons can 
be regarded as posing initial conditions for Starobinsky solutions of the anomaly 
driven deSitter expansion, which are thus dynamically eliminated by infrared effects 
of quantum gravity. The resulting landscape of hot universes treated within the 
cosmological bootstrap (the self-consistent back reaction of quantum matter) turns 
out to be limited to a bounded range of the cosmological constant, which rules out a 
well-known infrared catastrophe of the vanishing cosmological constant and suggests 
an ultimate solution to the problem of unboundedness of the cosmological action in 
Euclidean quantum gravity. 
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1. Introduction 

The ideas of quantum cosmology [U [2] and Euclidean quantum gravity [3], H] are again 
attracting attention. One of the reasons is the fact that the landscape of string vacua is 
too big [5] to hope that a reasonable selection mechanism can be successfully worked out 
within string theory itself. Thus, it is expected that other methods have to be invoked, 
at least some of them appealing to the construction of the cosmological wavefunction 
0, El E]. This quantum state arises as a result of quantum tunneling from the 
classically forbidden state of the gravitational field. The Hartle-Hawking wave function 
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of the Universe [3J H] describes nucleation of the de Sitter Universe from the Euclidean 
4-dimensional hemisphere, \I/hh ~ exp(— 5*e) = exp(37r/2GA) and has a negative action 
which diverges to — oo for the cosmological constant A — > 0. This implies a well-known 
infrared catastrophe of small cosmological constant — a vanishing A is infinitely more 
probable than any positive one. Alternative tunneling proposals for the wave function 
of the universe in the form of Linde [TO] or Vilenkin pj] give preference to big values of 
A, which opens the possibility for conclusions opposite to the Hartle-Hawking case. In 
particular, the inclusion of one-loop effects allows one to shift most probable values of the 
effective cosmological constant from zero to a narrow highly peaked range compatible 
with the requirements of inflation [12] . 

In this work we study the Hartle-Hawking prescription of the Euclidean path 
integration taking into account essentially nonlocal quantum effects mediated by 
nonlocal terms of non-vacuum nature. The core of our suggestion is a simple observation 
that the presence of radiation implies a statistical ensemble described by the density 
matrix, rather than a pure state assumed in [TJ [8] . Density matrix in Euclidean quantum 
gravity [13J, p[ip,ip'], originates from an instanton with two disjoint boundaries E and 
E' associated respectively with its two entries, see FigHJ Note that mixed nature of 
the quantum state is fundamental and the origin of impurity is not caused by coarse 
graining or tracing out environmental degrees of freedom. 
In contrast, the pure density matrix of the Hartle-Hawking state corresponds to the 



E E' 




Figure 1. Picture of instanton representing the density matrix. Dashed lines 
depict the Lorentzian Universe nucleating from the instanton at the minimal 
surfaces E and 

situation when the instanton bridge between E and E' is broken, so that topologically the 
instanton is a union of two disjoint hemispheres. Each of the half-instantons smoothly 
closes up at its pole which is a regular internal point of the Euclidean spacetime ball, 
see FigfS] — a picture illustrating the factorization of p = |^hh)(^hh|- 

When calculated in the saddle-point approximation the density matrix automat- 
ically gives rise to radiation whose thermal fluctuations destroy the Hartle-Hawking 
instanton. Namely, the radiation stress tensor prevents the half-instantons of the above 
type from closing and, thus, forces the tubular structure on the full instanton support- 
ing the thermodynamical nature of the physical state. The existence of radiation, in 
its turn, naturally follows from the partition function of this state. The partition func- 
tion originates from integrating out the field p in the coincidence limit ip' = tp. (This 
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Figure 2. Density matrix of the pure Hartle-Hawking state represented by the 
union of two vacuum instantons. 



procedure is equivalent to taking the trace of the operator exp(—j3H) inherent in the 
calculation of the partition function for a system with the Hamiltonian H.) This cor- 
responds to the identification of E' and E, so that the underlying instanton acquires 
toroidal topology. Its points are labeled by the periodically identified Euclidean time, a 
period being related to the inverse temperature of the quasi-equilibrium radiation. The 
back reaction of this radiation supports the instanton geometry in which this radiation 
exists, and we derive the equation which makes this bootstrap consistent. 

We show that for the radiation of conformally-invariant fields the analytical and 
numerical solution of the bootstrap equations yields a set of instantons — a landscape 
— only in the bounded range of A, 



This set consists of the countable sequence of one-parameter families of instantons which 
we call garlands, labeled by the number k = 1,2,3, ... of their elementary links. Each 
of these families spans a continuous subset, A^ n < A < AmL, belonging to ([Tj). These 
subsets of monotonically decreasing lengths Ami x — A^ m ~ l/k 4 do not overlap and their 
sequence has an accumulation point at the upper boundary A max of the range ([I]). Each 
of the instanton families at its upper boundary Aml x saturates with the static Einstein 
Universe filled by a hot equilibrium radiation with the temperature Tqa ~ mpj In k 2 , 
k ^> 1, and having the negative decreasing with k action Iq ~ — In 3 k 2 /k 2 . Remarkably, 
all values of A below the range (OQ), A < A min , are completely ruled out either because 
of the absence of instanton solutions or because of their infinitely large positive action. 

The details of our approach have recently been presented in [TJ], while here we 
discuss its basic points and, in particular, establish its relation to the well-known 
Starobinsky model [15] of the self-consistent deSitter expansion driven by the conformal 
anomaly of quantum fields. 

2. The effective action, generalized Priedmann equations and bootstrap 

As shown in [2] the effective action of the cosmological model with a generic spatially 



closed FRW metric ds 2 = N 2 (r)dr 2 + a 2 (r)d 2 n {3) = a 2 (dr] 2 + d 2 n^) sourced by 



A 



mm 



< A < A 



max' 



(1) 



conformal quantum matter has the form 
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+ B J + drN/a + f(2 J + drN/a^j , (2) 

where a(r) is the cosmological radius, N(t) is the lapse function, H 2 = A/3 and 
integration runs between two turning points at t±, a(r±) = 0. Here the first line is 
the classical part, the second line is the contribution of the conformal transformation 
to the metric of the static instanton ds 2 = drf + d 2 Vt^ (77 is the conformal time) 
and the last line is the one-loop action on this static instanton. The conformal 
contribution ri_i oop [g] — ri_i oop [g] is determined by the coefficients of QR, the Gauss- 
Bonnet invariant E = R 2 —4R 2 +R 2 and Weyl tensor term in the conformal anomaly 



g/j.uSTi-ioop/Sg^ = g 1 / 2 (anR + (5E + ^C 2 $) /4(47r) 2 . Specifically this contribution can 



be obtained by the technique of [16J; it contains higher- derivative terms ~ d 2 which 
produce ghost instabilities in solutions of effective equations. However, such terms are 
proportional to the coefficient a which can be put to zero by adding the following finite 
local counterterm 

T R [g] = rV loop [<?] + — J d'xg'l 2 R\g). (3) 

This ghost-avoidance renormalization is justified by the requirement of consistency of 
the theory at the quantum level. The contribution T R [g] — I\r[<?] to the renormalized 
action then gives the second line of ([2]) with B = 3/3/4. 

The static instanton with a period 770 playing the role of inverse temperature 
contributes ri_i oop [p] = E r] + F(r) ), where the vacuum energy E and free energy 
F(j]o) are the typical boson and fermion sums over field oscillators with energies wona 
unit 3-sphere E = i^^cj/2, F(r)o) = ±^ w ln(l =pe~ wr?0 ). The renormalization ([3]) 
which should be applied also to ri_i oop [p] modifies E , so that T R [g] = C r] + F(rj ), 
Cq = Eq + 3a/ 16. This gives the third line of Eq.(j2J) with Cq = B/2. This universal 
relation between Co and B = 3/3/4 follows from the known anomaly coefficients p2] and 
the Casimir energy in a static universe [18] for scalar, Weyl spinor and vector fields. 

The Euclidean Friedmann equation looks now as 

C = B/2 + F'fo), Vo = 2^ + dr/a(r). (5) 

The contribution of the nonlocal F(r)o) in (j2J) reduces here to the radiation constant 
C as a nonlocal functional of a{r), determined by the bootstrap equation (j3j), F'(j]o) = 
dF(i]o)/di]o > being the energy of a hot gas of particles, which adds to their vacuum 
energy B/2. 

Periodic instanton solutions of Eqs.(j4])-(l5]) exist only inside the curvilinear wedge of 
(if 2 , C)-plane between bold segments of the upper hyperbolic boundary and the lower 
straight line boundary of Figj3l 

4CH 2 <1, C>B-B 2 H 2 , BH 2 < 1/2. (6) 
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Figure 3. Instanton domain in the (H 2 , C)-plane. Garland families are shown 
for k = 1, 2, 3, 4. Their sequence accumulates at the critical point (1/25, B/2). 



Below this domain the solutions are either complex and aperiodic or suppressed by 
infinite positive Euclidean action. Indeed, a smooth Hartle-Hawking instanton with 
a_ = yields r/ — > oo in view of (jSJ), so that F(r} ) ~ F'iVo) — > 0. Therefore, its 
on-shell action 



due to 5 > diverges to +oo at a_ = and completely rules out pure-state instantons 
[T4"] . For the instanton garlands, obtained by glueing together into a torus fe copies of 
a simple instanton [19] the formalism is the same as above except that the conformal 
time in (jSJ) and the integral term of (jjTJ) should be multiplied by fc. As in the case 
of k — 1, garland families interpolate between the lower and upper boundaries of 
They exist for all k, 1 < k < oo, and their infinite sequence accumulates at the critical 
point C = B/2, H 2 = 1/25, with the on-shell action T ( Q k) ~ -5 In 3 k 2 /4k 2 n 2 which, in 
contrast to the tree-level garlands of [19] is not additive in k. 

3. Conclusions: Euclidean quantum gravity and the Starobinsky model 

Elimination of the Hartle-Hawking vacuum instantons implies also ruling out well-known 
solutions in the Starobinsky model of the anomaly (and vacuum energy) driven deSitter 
expansion [15] . Such solutions (generalized to the case of a nonzero "bare" cosmological 
constant 3H 2 ) can be obtained by the the usual Wick rotation from the solutions of (j4j) 
with the thermal contribution switched off. The corresponding Euclidean Friedmann 
equation reads 



and has a generic solution of the form a = sin hr/h, 1/ a 2 — a 2 / a 2 = h 2 , with the following 
two values of the Hubble parameter h = h±, h 2 ± = (1/B)(1 ± y/1 - 2BH 2 ). For H = 






(8) 
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this is exactly the Euclidean version of the Starobinsky's solution [TS] with h 2 + = 2/B. 
For larger H 2 < 1/2B we have two families of exactly deSitter instantons which can be 
regarded as initial conditions for the "generalized" solutions of Starobinsky. However, 
all of them are ruled out by their infinite positive Euclidean action. 

What remains is a quasi-thermal ensemble of non-vacuum models in the bounded 
cosmological constant range (JT]) with A min > and A max = 3/25 and with a finite 
value of their effective Euclidean action. This implies the elimination of the infrared 
catastrophe of A — > and the ultimate solution to the problem of unboundedness of 
the on-shell cosmological action in Euclidean quantum gravity. As a byproduct, this 
suggests also strong constraints on the cosmological constant apparently applicable to 
the resolution of the cosmological constant problem. 
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